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It is generally assumed that gas in a galactic disk follows closely non self- 
intersecting periodic stellar orbits. In order to test this common assumption, we have 
performed MHD simulations of a galactic-like disk under the influence of a spiral 
galactic potential. We also have calculated the actual orbit of a gas parcel and com- 
pared it to stable periodic stellar orbits in the same galactic potential and position. 
We found that the gaseous orbits approach periodic stellar orbits far from the major 
orbital resonances only. Gas orbits initialized at a given galactocentric distance but at 
different azimuths can be different, and scattering is conspicuous at certain galacto- 
centric radii. Also, in contrast to the stellar behaviour, near the 4:1 (or higher order) 
resonance the gas follows nearly circular orbits, with much shorter radial excursions 
than the stars. Also, since the gas does not settle into a steady state, the gaseous 
orbits do not necessarily close on themselves. 

Key words: Galaxy: disk - Galaxy: kinematics and dynamics - Galaxy: structure - 
galaxies: spiral structure - MHD 



1 INTRODUCTION 

Differences between stellar orbital dynamics and gas dynam- 
ics in disk galaxies play a fundamental role in the interpre- 
tation (or misinterpretation) of kinematics observations in 
order to deduce gener al physical properties of galaxies. One 
example (Gomez 2006) is that the inferred quantity of dark 
matter in galaxies can be affected by inconsistencies between 
the well-studied gas rotation curve and the relatively poorly 
studied stellar disc rotation curve. 

Several observational and theoretical studies have 
been devoted to the investigation of the relationship 
between the gas streamlines and the corresponding stel- 
lar orbits. In particular, it is customary to identify the 
non self-intersecting stable periodic orbits of a given 
astronomical system, with plausibl e regions for ^as 
strea ms to settle (e.g. iKalnaisI Il973l: Simonson Mader 
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From resonances emerge strong orbital families like 
the xl type that manage to survive in models up to ap- 
proximately the 4:1 resonance. After that periodic orbits 
from many families are found. This means that at reso- 
nances new families of periodic orbits start existing. At res- 
onances wejnaY_encou^^ than o n e family of periodic 
orbits dContoDOulos Grosb"0l Il986l : iHeller Shlosmanl 
Il996l : iKaufma nn Patsisll2005h . This is in part intrinsic to 
the problem, due to the non-linear character of the pertur- 
bations. In normal non-barred spiral galaxies with a moder- 
ate non-axisymmetric component, such as spiral arms, the 
stellar orbits and gaseous flow can be highly non-linear (this 
means, the response of a dynamical system to a perturbation 
depends non- linearly on the amplitude of the perturbation), 
even far from resonances. Also, phenomena like gas shocks or 
stellar chaos are f ound even in the case of low-amplitude spi- 
ral pe rturbations (|Pichardo et al.|[2003l : IPerez- Villegas et al.l 
I2OI2I I. 

Large scale galactic gaseous and stellar flows present 
various important differences in their kinematics. Some in- 
teresting known differences between the response of gas 
and stars to an imposed galactic potential model are due 
to the fact that gas is very responsive and suffers violent 
shocks, unlike stars. In resonant regions, it has been sug- 
gested that sho cks occur close to the intersections of stable 
periodic orbits (Schwarz 1981; van Albad a Sander si [198^ : 
lAthanassoula 1992b). Gas is attracted to the spiral arms 
by their gravitational field and, unlike stars, it responds to 
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pressure gradients. In non-barred galaxies with large bulges 
or thick disks (such as early spirals), stars are supported 
not only by rotation but by strong radial velocities or veloc- 
ity dispersion, while gas settles down in much thinner disks 
supported by rotation. The later the morphological type of 
the gal axy, the m ore similar are gas and stellar kinematics 
(Beck man, Zurita Beltran 2004). In this work we study 
the interplay between gas flow and stellar orbital dynami- 
cal properties of spiral galaxies. We have calculated a set of 
MHD and particle simulations in a static background grav- 
itationa l potent i al ba sed on a full 3D spiral galactic model 
(Pichardo et al.ll2003l ). This includes spiral arms different 
in nature to the classic tight winding approximation, in or- 
der to contribute to the understanding of the relationship 
between stellar orbits and gas orbits in non-barred spiral 
galaxies, and the relation with periodic orbits inside and 
outside resonance regions. 

This paper is organized as follows. In 321 the galactic 
potential used to compute the MHD and stellar orbits and 
the initial simulation set up is described. In ^ the stellar 
and MHD orbits obtained with the spiral potential are pre- 
sented and compared with each other. In 311 we present a 
brief summary and discussion of our results. 



2 NUMERICAL MODEL 

For the purpose of this work we have produced MHD simula- 
tions under the effect of a detailed model of the background 
galactic potential. We describe next the potential model and 
the performed simulations. 



2.1 The Galactic Background Model 

For the orbital study we have employed a detailed galac- 
tic semi-analytic model, resembling a Milky- Way- like poten- 
tial. Rather than using a simple ad hoc model for a three- 
dimensional spiral perturbation (such as a cosine function), 
we constructed a three-dimensional mass distribution for the 
spiral arms and derive their gravitational potential and force 
fields from previously known results in potential theory. The 
spiral arms distribution consists of inhomogeneous oblate 
spheroids superposed along a given spiral locus. It is phys- 
ically simple, with continuous derivatives and density laws. 
In principle, the dimensions and mass density of the oblate 
spheroids will depend on the type of spiral arms that are 
modeled, gaseous or stellar. The model, called PERLAS, is de- 
scribed in detail in IPichardo et al.l (|2003i ) . Comparisons of 
the model with other models have been already published, 
in particular for the most relevant parts of our stu dy, the 
spiral arms perturbation ( Antoja et al. 2009; Martos et al.l 
[20041 : IPichardo et al.ll2003l l . The corresponding observation- 
ally motivated parameters of the model used for this work 
are presented in Table [1] The spiral perturbation is a 3D 
steady two- a rmed model that traces the locus reported by 
[Drimmel Spergel ( 2001), using K-band observations. The 
solar radius, at 8.5 kpc, is close to the spiral arms 4:1 res- 
onance. The mean force ratio between the arms and the 
axisymmetric background is a round 10%, which is in agree- 
ment with the estimations by Patsis et al. (1991) for Milky 
Way type galaxies. At the solar circle, the radial compo- 
nent of the force due to the spiral arms is 4.4% of the mean 
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Figure 1. Density distribution after 1 Gyr (top), 2.5 Gyr (mid- 
dle), and 4 Gyr (bottom). Although the computational grid ex- 
tends to r = 22 kpc, only the inner region is shown. The dashed 
circles show the position of the inner Lindblad, 4:1, corotation, 
-4:1 and outer Lindblad resonances. Gas rotates in the clockwise 
direction (counterclockwise outside the corotation resonance). 
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Table 1. Non-axisymmetric galactic model (|Picliardo et al.ir2003l ) 



Parameter 



Value 



Reference 



Spiral arms locus 

Spiral arms pitch angle 

Spiral arms external limit 

Spiral arms: exp. with scale-length 

Spiral arms force contrast 

Spiral arms pattern speed {VLp) 



Bi-symmetric (Logthm) 
15.5° 
12kpc 
2.5 kpc 
- 10% 
20 km s~^ kpc 



Churchwell et al. r2Q09 
Drimmel (2000) 
Drimmel (2000) 
Disk based 
Vot sis et al . (1991) 
Martos et al. (^2004) 



axisymmetric background. The self-consistency of the spiral 
arms has been tested through the reinfor cement of the spi- 
ral potential by the s table periodic orbits |Patsis et al.ll 199 ll : 
IPichardo et~ar]l2003l ). 

2.2 The MHD Setup 

The initial setup of the MHD simulations consists of 
a gaseous disk with a density profile given by n(r) = 
noexp[— (r — ro)/rd], where no — 1.1 cm~^, ro — 8 kpc, 
and Td — 15 kpc. The gas behaves isothermally with a 
temperature T = 8000 °K and is permeated by a mag- 
netic field, initially in the azimuthal direction, with inten- 
sity B{r) = 5oexp[— (r — ro)/rB] where Bq = 5 /xG and 
tb — 25 kpc. The disk is initially in rotational equilibrium 
between the centrifugal force, the thermal and magnetic 
pressures, magnetic tensi on and a background axisymmet- 
ric galactic potential (All en San tillan 1991) consisting of 
a bulge, stellar disk and halo. 

This equilibrium is perturbed by the PERLAS spiral arm 
potential, which rotates with a pattern speed of Qp = 
20 km s~^ kpc~^. This pattern speed and background po- 
tential place the inner Lindblad, 4:1, corotation, -4:1, and 
outer Lindblad resonances at 2.81, 6.98, 10.92, 14.37 and 
17.72 kpc, respectively. 

We solved the MH D equations with the ZEUS code 
(jStone Normanlll992al lbll. which is a finite difference, time 
explicit, operator split, Eulerian code for ideal magnetohy- 
drodynamics. We used a 2D grid in cylindrical geometry, 
with R G [1.5 kpc, 22 kpc] and a full circle in the azimuthal 
coordinate, 0, using 750 x 1500 grid points. Both boundary 
conditions in the radial direction were outflowing. All the 
calculations are performed in the spiral pattern reference 
frame. No self- gravity of the gas was considered. 

The simulation starts with the gas in circular orbits in 
equilibrium with the background galactic potential, thermal 
and magnetic pressure gradients and magnetic tension. After 
the perturbation is activated, the gas very rapidly settles 
into a spiral pattern, with two arms (plus an oval-shaped 
rin£l) inside the inner Lindblad resonance (ILR) and four 
arms outside, up to the 4:1 resonance (see fig. [I]). The inner 
pair of arms follows a ^ 9° pitch angle, while the two pairs 
outside the ILR follow a ^ 9° and ~ 13° pitch angles; in 
contrast, the perturbing arm potential follows a 15.5° pitch. 
The exact position and pitch angle of the gaseous spiral arms 
oscillate slightly around the above quoted values. 

^ This ring at r < 2 kpc falls in a region where the dynamics 
should be dominated by a galactic bar, which is not included in 
this work. Therefore, it is disregarded. 
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Figure 2. Comparison of the stellar (contours) and gaseous den- 
sities (grayscale) in the simulation. The stellar density is that 
corresponding to the background axisymmetric galactic disk plus 
the (imposed) spiral arms potential, while the gas density corre- 
sponds to to that of a snapshot taken 1 Gyr after the start of the 
simulation. Since the model corresponds to a trailing spiral, the 
gas flows down from the top of the plot. The differences in number 
and position of the stellar and gaseous arms are apparent. 

It is noticeable the presence of a n MHD instability at 
corotation radius (further described in lMartosl l2013) , which 
begins to develop 2.8 Gyr after the start of the simulation. 
This instability is present in a variety of simulations with 
different magnetic field intensities, and absent in purely hy- 
drodynamical ones. With exception of this instability, ab- 
sence of a magnetic field does not affect our conclusions. 

Figure [2] compares the imposed stellar density with the 
response gas density at a time 1 Gyr after the start of the 
simulation. The stellar density is obtained by taking the 
Laplacian of the potentials corresponding to the axisymmet- 
ric background disk and spiral arm distribution, as given by 
the PERLAS model. In this figure, the differences in position 
and pitch of stellar and gaseous arms are apparent. 



3 ORBITS IN THE SPIRAL POTENTIAL 

Figure [3] shows a comparison of stable periodic stellar 
and gaseous orbits in the background axisymmetric+spiral 
galactic potential. That figure also shows the positions of the 
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Figure 3. Comparison of some stellar {red) and gaseous {green) 
orbits near major resonances. The location of the inner Lindblad 
(ILR), 4:1 and corotation are indicated by dashed circles. The 
gaseous orbit near the ILR starts similar to the stellar orbit, but 
it rapidly decays to a ring. Around other resonances, the gas 
follows a near circular orbit, even as the stellar orbits experience 
large radial excursions. As reference, the gas density {grayscale) 
and the locus of the stellar arms {black line) are also shown. 



major resonances and the locus of the spiral arms as defined 
in the PERLAS model. Since ZEUS is an Eulerian code, it does 
not provide the actual path a gas parcel follows. Therefore, 
in order to reconstruct a gaseous orbit, we integrated the 
velocity field interpolating in space and time between data 
files for an array of initial positions, following the trajecto- 
ries of individual gas parcels in the resulting velocity field. 
The integration was performed over 4 Gyr, starting 1 Gyr 
into the simulation in order to avoid the initial settling of 
the gas into the spiral pattern. If during that time the gas 
parcel reaches either of the numerical boundaries (at = 1.5 
and 22kpc), we stop the integration. The stellar orbits were 
obtained by direct integration in the potential, locating the 
periodic ones using a Newton-Raphson algorithm. 

When comparing the stellar and gaseous orbits pre- 
sented in Figure O it is apparent that some of them are 
similar, but also some may be quite different, even when 
both gas and stars are subjected to the same gravitational 
potential. 

Stars and gas follow orbits that are similar in the regions 
between the resonances. In Figure [4l we show some of the 
orbits calculated in similar regions; we also show the stellar 
mass and gas densities for comparison with the orbit. In the 
case of the stellar orbits, it is noticeable how the familiar 
oval shape (associated with the ILR) gradually morphs into 
a rounded square as the orbit approaches the 4:1 resonance. 
These last orbits are less supportive of the spiral pattern, in 
accordance with analytical theory. 

Although the gas orbits are complex and quite sensi- 
tive to the initial position, they follow a similar evolution to 
the stars when regions away from the ILR are considered. 
It is noticeable that the gaseous orbits develop features that 



give rise to the second pair of spiral arms (apparent in the 
rightmost column of fig.[4j). As opposed to the stars, the gas 
may change direction rapidly by developing oblique shocks, 
which generate pressure gradients that further change the 
path of the orbit. Between the ILR and 4:1 resonance, the 
stellar and gaseous orbits are similar in the radial range they 
span, but the shape differs when the second pair of gaseous 
arms becomes stronger, since the shocks at the arm posi- 
tions significantly change the path of the gas parcel. Just 
outside the 4:1 resonance, at the same Jacobi constant, we 
have more than one main family of periodic orbits coexist- 
ing, some of them stable, but the most of them unstable 
(see fig. 7 in IContopoulos Grosb^ Il986l ). Most of stars 
approach an orbit with a four- fold symmetry; the gas, on 
the other hand, selects a twofold symmetric orbit related to 
the 2:1 coexisting periodic orbits family, associated with the 
only pair of arms found in the gas in this region. 

Near the resonances, the gas and stars follow quite dif- 
ferent orbits (fig. El). While the stars follow the familiar geo- 
metrical shapes, the gas is unable to follow the large depar- 
tures from circularity of the stellar orbits. While the stars 
follow oval-shaped periodic orbits near the ILR, the gas ac- 
tually tries to avoid that region. Integrations starting near 
or inside the ILR decay to the ring/spur structure just inside 
2 kpc. This structure is not expected to appear in real galax- 
ies, since usually it is affected by the starting position of the 
spiral arms and should be strongly influenced by a galactic 
bar not included in the present model). In fact, in order to 
find a gaseous orbit that resembles those found in the stars 
near the ILR, we must look for an orbit that remains out- 
side the ILR region during its time evolution (the top row in 
fig.m for example). Also, these gaseous orbits present some 
dispersion, since they are not quite periodic and actually 
cross themselves. This is possible since the simulation does 
not really reach a steady state, but instead settles on a pe- 
riodic cycle, with the arms oscillating slightly around given 
positions. 

Near the 4:1 resonance, the periodic stellar orbit follows 
the familiar rounded-square, while the gas actually remains 
close to a circular orbit, with deviations from circularity of 
less than 400 pc (although these deviations are systematic 
and allow for one pair of the gaseous spiral arms to extend 
beyond this resonance, as opposed to the stellar pattern). 
The strong evolution of the gaseous orbits around the ILR 
is absent near the 4:1 resonance. 

Further out, the expected banana-shaped orbits appear 
at some positions, although the gaseous orbits are greatly 
distorted due to the turbulence related to the instability at 
corotation (fig. [6]). The strong velocity perturbations related 
to this instability appear to generate points that act as at- 
tractors. Future work will try to determine whether these 
positions are suitable for star formation. Outside corota- 
tion, the gas orbits are nearly circular again, with deviations 
less than 0.5 kpc from the initial radius. This is not surpris- 
ing since the underlying spiral perturbation is already weak, 
representing only 1% of the radial axisymmetric force at the 
-4:1 resona nce (the arm potent ial tapers-off at a radius of 
12 kpc, see lPichardo et al-lliooi ). 
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Figure 4. Comparison of stellar and gaseous orbits between resonances. The first two columns show stellar orbits in cartesian {first 
column) and polar coordinates {second column) plotted over the stellar mass density {grayscale) corresponding to the PERLAS potential 
model. The last two columns show gaseous orbits, also in cartesian and polar coordinates, plotted over the gas density {grayscale) 
averaged over the 4 Gyr period of the orbit integration. In all cases, the initial position of the integration is marked by the small open 
circle. Since the model corresponds to a trailing spiral, the sense of the gaseous orbits is clockwise in the cartesian coordinate plots, and 
down from the top in the polar coordinate ones; only stellar orbits with the same sense of rotation were considered. The position of the 
inner Lindblad, 4:1 and corotation resonances is shown {dashed circles). 



4 DISCUSSION AND CONCLUSIONS 

We performed MHD simulations of a gaseous disk subjected 
to the PERLAS spiral potential (disregarding the effect of a 
galactic bar, which will be explored in the future). We then 
studied the actual path the gas follows as it rotates around 
the galaxy in order to compare it with the stable periodic 
stellar orbits existing in that same galactic potential in order 
to test the frequently stated assumption that the gas should 



follow orbits close to the periodic stellar orbits that do not 
cross themselves. We found that the gas does not always 
behave that way. 

The most obvious difference in the stellar and gas be- 
haviour is that the gas res ponds to the two-arm poten- 
tial with four spir al arms Jshu. Milione k RobertsI Il973l : 
iMartos et al.l lioM ), organized in two pairs, each with a 
tighter pitch angle than the underlying potential. The four 
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Figure 5. Similar to fig. [U for orbits near the resonances. 
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Figure 6. Gaseous orbits near the corotation resonance. 



gaseous arms are well defined between the ILR and the 
4:1 resonance, more or less straddling the stellar arms in 
similar fashion to the se c ondar y compressions described by 
IShu. Milione Robertsl (|l973l ) (see fig. [2|). Outside the 4:1 
resonance, the gaseous arms follow the stellar arms more 
closely. The gas response to stellar arms depend on the de- 
tails of the imposed potential (compare, for example, with 



IPatsis et al.lll997l and lVorobvovll2006l l and will be explored 
in future work. 

As opposed to stars, the gas may develop shocks at 
the arm positions, which allow its orbits (understood as the 
path a gas parcel follows) to change direction abruptly. This 
behaviour generates cusps in the gaseous orbits that are not 
present in the stellar orbits. 
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Setting aside the extra cusps in the gaseous orbits, the 
stellar and gaseous orbits are most similar in between the 
resonances, when the general shape and radial range are 
considered. These are the regions where the periodic stellar 
orbits are rounder and the radial excursions are smaller, al- 
lowing the gas to more closely follow the stellar orbits since 
it will be less influenced by forces other than gravity. Near 
the resonances, however, the stars suffer large radial excur- 
sions that imply regions with small curvature radii so that 
the orbit may close on itself (see fig.[5l). A gas parcel cannot 
follow such an orbit0 and either looses angular momentum 
and moves away from the resonance (such as the orbit near 
the ILR), or settles in an orbit with only small radial excur- 
sions. The gravitational forcing that maintains the resonance 
is still present; so, even as the gas follows a nearly circular 
orbit, its velocity is not uniform, thus generating pressure 
waves that allow the spiral arm to extend across the res- 
onant radius; this is apparent, for example, across the 4:1 
resonance. 

Although banana-type orbits should be present at the 
corotation radius, the gas only develops such orbits just be- 
yond that resonance. Generally speaking, the gas experi- 
ences changes always in the outer side of the resonant radii, 
for example, the second pair of arms start and end just out- 
side the ILR and 4:1 resonance, respectively. This particular 
element will be further explored in future work. 

When stating that the gas should follow stable periodic 
stellar orbits that do not cross, a key assumption is that the 
gas falls into a steady state. Small oscillations of the gaseous 
arms (an effect exacerbated when t he vertical structure 
of the gaseous disk is considered, see iGomez k Coxll2004l ) 
and magnetohydrodynamic instabilities (like that present at 
corotation, see Martos 2013), distort the gaseous orbits so 
that they may cross themselves. Other phenomena not con- 
sidered here, like star formation and feedback, will further 
disrupt the path a gas parcel should follow. 

On the other hand, in chaotic regions, such as the coro- 
tation zone, the phase space for stable periodic orbits is 
severely reduced. Gas dynamics behaviour in this region can 
not be attributable to stable sim ple paths given by periodic 
orbits (|ChatzoDOulos et"al]|201lh . In these regions, chaotic 
stellar orbits might be the ones supporting large scale struc- 
tures, and gas might not be able to follow chaotic orbits 
since their trajectories are, in general, self-intersecting. 

Since observations of the spiral structure of the Milky 
Way are well fitted by four arms in the gas and two in the 
stars, the orbits that support such structures cannot be the 
same. The differences in the orbital behaviour of gas and 
stars shown here highlights the difficulty of extrapolating the 
results of stellar dynamics to gas dynamics. Even if the mo- 
tion is dominated by gravity, other physical processes may 
have a strong impact on the overall behaviour of a galactic 
disk. 



^ We have also found that there is a limit to the ellipticity of a 2:1 
orbit (oval ones) at approximately 0.4, above which gas cannot 
settle down readily on the orbit because of the presence of strong 
shocks at the apocenter. 
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